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I. INTRODUCTION 



A flexible polymer chain in a dilute solution is influenced by both hydrophobic interac- 
tions between the monomers and the excluded volume effect. The attractive interactions are 
neglected at high temperatures or in a good solvent, but become significant as the temper- 
ature T is lowered. As T reaches a special temperature 9, the linear polymer undergoes an 
abrupt change from an expanded conformation for T > 9 to a fully compact conformation 
for T < Long polymer in a good solvent is a critical system, and the collapse transition 
at T = 9 has been identified as a tricritical transition 3 - 1 ^. The 9 point behavior is well- 
described by self-avoiding walks with attractive interaction energy assigned for each pair of 
nonbonded nearest-neighbor monomers. The tricritical exponents take the mean-field values 
for d > 3, and there are logarithmic corrections at d = 3^-. A great deal of studies have 
been performed to understand the nature of the collapse transition in two dimensions^"—, 
which is expected to exhibit much more non-trivial behavior than its higher dimensional 
counterparts. 

In this work, we study the collapse transition of a polymer on a square lattice, with both 
nearest-neighbor (NN) and next nearest-neighbor (NNN) interactions present, by calculat- 
ing the exact partition functions up to chain length N = 36. We estimate the crossover 
exponent and the transition temperature from the zeros of the partition function, and also 
from the specific heat. Although the method of partition function zeros became one of the 
most popular tools for studying the critical phenomena with the advancement of computa- 
tional power— there are few works where partition function zeros of lattice polymers were 
calculated. For examples, exact partition function zeros were computed for the simple-cubic 
lattice up to chain length 13 36 , for the face-centered lattice up to chain length 9 37 , and for 
the square lattice up to chain length 36 28,29 . Only NN interactions were present in these 
works. In fact, the current work is the first instance where a square- lattice polymer with 
NNN interactions is ever studied. It was only on a hexagonal lattice that models with NNN 
interactions were studied previously -"— . 

By introducing NNN interactions, the transition behavior is much more pronounced than 
that of the model only with NN interactions 29 . The results suggest that the model belongs 
to the same universality class as the one described by the model with only NN interactions. 
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II. THE NUMBER OF CONFORMATIONS 



Conformations of a polymer chain with N monomers are modeled as a two-dimensional 
self-avoiding chain of length N on a square lattice. The position of the monomer i is given 
by rj = (k,l), where integers k and I are the Cartesian coordinates relative to an arbitrary 
origin. Chain connectivity requires |r, — rj+i| = 1, i.e., bond length is unity. Due to the 
excluded volume, there can be no more than one monomer on each lattice site, ^ Tj 
for % ^ j. The attractive hydrophobic interaction is incorporated by assigning the energies 
—ei < and — e 2 < for each non-bonded NN and NNN contact between monomers. The 
resulting Hamiltonian is 

H = -e 1 J2 A(ri, r,) - e 2 ]T A(r<, r,), (1) 



i<j i<j 



where 



and 



.1 if \i — 7 1 > 1 and |r,- — r,-| = 1, 
A(r 4 ,r,)=<( ' Jl lf Jl (2) 

otherwise, 

,1 if |r,- — r,-| = v2, 
A(r,r,)H ' l 11 V (3) 

otherwise. 

The result when only NN interactions are present, corresponding to the 9 point^, can be 
reproduced by putting e 2 = 0. We consider the case with e\ = e 2 = e. The energy of the 
system is then E = —e{K\ + K2) = — eK, where K\ and K2 are the number of contacts 
between NN and NNN monomers, respectively. 

Here we define the reduced number of conformations ojn(K), where conformations related 
by rigid rotations, reflections, and translations are regarded as equivalent, and counted only 
once. On the other hand, due to an assumption that the polymer chain has an intrinsic 
direction, the conformations with reverse labels i •<-)■ N — % + 1 for all (i = 1, 2, • • • , N) are 
considered distinct. It is easy to see that the total number of conformations generated by 
rotations and reflections from a given conformation is eight, except for the straight chain 
where the total number of conformations generated by rotations and reflections is four due 
to invariance with respect to reflection perpendicular to the chain. The total number of 
conformations Qn(K) is obtained from u>n(K) as follows: 

iu N (K)-A if if = 0, 

n N (K) = { (4) 

1u)n(K) otherwise. 



Thus, one can achieve about eight-fold reduction in the computing time by enumerating the 
reduced number of conformations un(K) instead of Qn(K)2&. We obtained un(K) up to 
A" = 36 by the help of a parallel algorithm classifying conformations by sizes of rectangles 
they span2£. 

III. PARTITION FUNCTION ZEROS IN THE COMPLEX 
TEMPERATURE PLANE 

Yang and Lee^ first introduced the concept of the partition function zeros in the complex 
fugacity plane, and found a mechanism for the occurrence of phase transitions in thermo- 
dynamic limit. Later, Fisher— showed that the partition function zeros in the complex 
temperature plane are very important in understanding phase transitions. For system ex- 
hibiting the temperature-driven phase transition, the locus of Fisher zeros forms a line and 
crosses the positive real axis in thermodynamic limit. The intersection point of the locus 
with the positive real axis corresponds to the critical temperature. The zeros closest to 
the positive real axis are called the first zeros, which approach the positive real axis as the 
system size increases. 

The partition function of our model is 



where y = exp(/3e) and (3 = 1/ksT. We see that since K is bounded, the partition function 
(JSJ) is a n-th order polynomial of y where n is the maximum value of K. The partition 
function zeros yi [i = 1,2, ••• ,n) are then obtained by solving the polynomial equation 
Z(y) = 0. The solution was found with mathematica. As can be seen from Fig. [TJ the 
first zeros approach the positive real axis in the complex temperature plane as polymer 
length increases. 

IV. THE SCALING BEHAVIOR AND THE CRITICAL EXPONENT 

Near the critical temperature T c , the radius of gyration (or the end-to-end distance) Rn 
of a polymer chain with A" monomers is generally expressed by the scaling theory^, 




(5) 



K 



(R%) ~ N 2v f{rN% 



(6) 
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where the reduced temperature is defined as r = |T — T c \ /T c and the scaling function f(x) 
behaves as follows: 



x (6/(d+2)-2v)/<p if x ^ 00) 

const. if x — >■ 0, (7) 

x (2/dr-to)/4> if X ^ -OO. 



The exponent z/ represents the geometrical properties of a polymer, and the crossover expo- 
nent describes how rapidly the system undergoes the transition as T approaches T c . The 
crossover exponent also describes how rapidly the first zeros approach the positive real 
axis as N increases^, 

lm[ yi (N)} ~ j\T* (8) 

where yi(N) is a first zero for a polymer chain with N monomers. In finite-size systems 
with even N, the crossover exponent is approximated as 

m _ ln{Im[y 1 (iV + 2)]/Im[y 1 (iV)]} 
W} ~ \n{{N + 2)/N} ' [) 

which reduces to the exact value of in N — > oo limit, estimated by using the Bulirsch-Stoer 
(BST) extrapolation^. We obtain 0.4422(14) for the crossover exponent as shown in Fig. [2, 
where the estimated error could be further reduced by removing unreliable data obtained 
from iV < 18. The error is estimated by examining the robustness of the extrapolated 
value with respect to perturbations of the data points, but it is not a statistically rigorous 
confidence level^ 1 ^. Therefore, we estimated the error by slightly changing the ratio of 
NNN and NN interactions, R = £2/^1, which we set to 1 in the current work. We change 
R by 0.5, and get = 0.428 for both R = 0.5 and 1.5. If we assume that R is irrelevant 
and combine the results for R = 0.5, 1.0, and 1.5, the resulting range of the crossover 
exponent is 0.428 < < 0.442. The result is consistent with the conjectured exact value of 
= 3/7 = 0.4286 obtained from hexagonal lattice with random annealed forbidden faces^, 
as well as our previous estimate from the model with NN interactions only, <fi = 0.422(12), 
suggesting that they belong to the same universality class. More extensive analysis for 
various values of R is postponed for a future study. 

Without additional information, we assumed the leading finite size correction to <fi is 
of order 0(N^ 1 ) when performing the BST procedure. We estimated the range of also 
by changing the leading exponent of the extraopolating function. With R = 1 fixed, we 



performed BST extrapolation with the leading finite size correction of order 0(iV -aJ ) with 
u = 0.5 and 1.5. We get <fi = 0.418 and 0.458 for u = 0.5 and 1.5 respectively, and combining 
these results with that for uj — 1.0, we get 0.418 < <p < 0.458, again consistent with both 
the conjectured exact value and the estimate from the model with NN interactions only. 
Again, there is no evidence that our model belongs to a universality class different from that 
of the model with NN interactions only. 

The real parts of the first zeros can be used to estimate the critical temperature y c , by 
estimating the point they approach in the limit of iV — > oo, 

Re[ yi (iV)]-y c ~JV-* (10) 

with the value of obtained above. The value of y c , obtained by extrapolating the data 
for even N with N > 18, is 1.3279(41), which corresponds to Tje = 3.526(39) (Fig. [3J. 
It is also shown in Fig. [1] along with the result for the model where only NN interactions 
are present^, corresponding to y c = 2.16(18) (T c = 1.30(17)). The transition temperature 
becomes much higher when additional attractive NNN interactions are included, which is to 
be expected. We obtain y c = 1.3288(41) with the conjectured exact value = 3/7 8 , which is 
not much different from the result above. As can be seen from Fig. [I], the transition behavior 
is much more visible when we introduce NNN interactions. 



V. SPECIFIC HEAT 



Now we estimate the critical temperature y c again by analyzing the behavior of the specific 
heat per monomer, for comparison with the result obtained from the partition function zeros. 
The specific heat per monomer is 
C(T,N) 1 dE 



e 2 N 



e 2 NdT 
/3 2 d 2 \nZ 

e 2 N dp 2 



(11) 



N 



Y, K K 2 n N {K)y 
J2 K n N (K)yK 



K 



Y.K^N{K)y 



K 



which is plotted in Fig. H] as a function of y for several values of N . The finite N approxima- 
tion of the transition point, y c (N), is obtained from the condition ^ = 0. We observe a peak 
around y ~ 1.5, which becomes sharper as N increases. By applying the BST extrapolation 
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to the finite-size scaling 

Vc{N) — y c (oo) ~ N~^, (12) 

we obtain the transition point y c (oo) = 1.265(19), equivalent to T c /e = 4.25(29), where the 
data for even N with 18 < iV < 36 were used. y c (N) is displayed in Fig. as a function 
of l/N®, along with the extrapolated value y c (oo). The current result is not drastically 
different from that obtained by the partition function zeros, but the precision is lower due 
to the fact that the specific heat is riddled by noisy contributions from zeros other than the 
first ones22. 

VI. DISCUSSIONS 

In this work, we studied the collapse transition of a square-lattice polymer with both 
NN and NNN interactions, by calculating the exact partition function zeros up to chain 
length N = 36. The crossover exponent <fi and the transition temperature T c were obtained 
by examining their scaling behavior with increasing chain length. We estimated T c also by 
calculating the specific heat from the exact partition function. Our results suggest that the 
polymer with both NN and NNN interactions on a square lattice belongs to the 9 universality 
class described by the model where only NN interactions are present, but by introducing 
NNN interactions, the transition behavior becomes more pronounced than the model with 
only NN interactions^. 
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FIG. 1. Positions of the first zeros in the first quadrant of the complex temperature (y = e^ e ) plane 
for TV = 10, 12, • • • , 36. Open circles indicate the results when both NN and NNN interactions are 
present, and open squares are those for the model with NN interactions only. Two dots indicated 
by arrows are the corresponding values of y c . 
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FIG. 2. The finite size approximations of the crossover exponent, (j>{N), are shown as a function 
of 1/N for even N with 10 < N < 18 (open circles) and N > 18 (solid circles). The value 
of (j) = 0.4422(14) for N — > oo (the open circle with an error bar) is estimated by the BST 
extrapolation for N > 18. 



11 



1.4 



1.3 



1.1 



0.9 



1 1 


1 1 








o 


— 


o 




o 


1 1 


o 

1 1 



0.0 



0.1 0.2 
1/A^ 



0.3 



0.4 



FIG. 3. The real parts of the first zeros are shown as a function of \/N^ for even N with 10 < 
N < 18 (open circles) and N > 18 (solid circles). The value of y c = 1.3279(41) (the open circle 
with an error bar) for N — > oo is estimated by the BST extrapolation for N > 18 with 4> = 0.4422. 
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FIG. 4. The specific heat for N = 20, 28, and 36 from bottom to top. 
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FIG. 5. The finite size approximation of y c obtained from the specific heat, y c (N), are shown as a 
function of for even N with 10 < N < 18 (open circles) and N > 18 (solid circles). The value 
of y c (oo) = 1.265(19) (the open circle with an error bar) is estimated by the BST extrapolation for 
N > 18 with <p = 0.4422. 
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